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We discuss the available data for the differential and the total cross section for the photodisinte- 
gration of 3 He and 3 H and the corresponding inverse reactions below E~, — 100 MeV by comparing 
with our calculations using realistic NN interactions. The theoretical results agree within the er- 
rorbars with the data for the total cross sections. Excellent agreement is achieved for the angular 
distribution in case of 3 He, whereas for 3 H a discrepancy between theory and experiment is found. 
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I. INTRODUCTION 

Over the last decades the photodisintegration of 3 He 
and 3 H and the corresponding inverse reactions have been 
investigated experimentally and theoretically with con- 
siderable interest. There have been a lot of experiments 
usingdiffercnt techni ques for the photodisintegration of 
3 He [p]^l| and 3 H |22}-|2"§|l or the inverse reaction, respec- 
tively. Despite the many investigations, there are incon- 
sistencies between the data up to 30% in the magnitudes 
of the cross sections. 

Early theoretical calculations were restricted to phe- 
nomenological interactions and various approximations 
in the bound state wave function and the scattering states 
(for a list of References see |2?J ) . The first consistent cal- 
culation for both the initial and the final state was done 
by Gibson and Lehman j30). They solved Faddeev-type 
Alt-Grassberger-Sandhas (AGS) equations |3l]] using Ya- 
maguchi interactions and taking into account only the El 
contributions of the electromagnetic interaction. 

Attempts to use realistic interactions are the ones by 
Aufleger and Drechsel |32| and by Craver et al. j33j. In 
both calculations higher multipoles were considered, but 
the three-body scattering state was not treated exactly. 
The unusual energy dependence of the cross section pos- 
tulated by Craver et al. has never been confirmed by 
any other calculation. Klepacki et al. [p9| also used a 
realistic interaction, however, in plane wave impulse ap- 
proximation. King et al. Ea ] performed an effective two- 
body direct capture calculation with the initial state be- 
ing treated as a plane wave, or as a scattering state gen- 
erated from an optical potential. 

The very-low-energy n-d radiative capture process 
is dominated by the magnetic dipole (Ml) transition, 
and has been studied by several authors |34 35 in 
configuration-space with inclusion of three-body forces, 
final state interaction (FSI), and explicit meson ex- 
change currents (MEC) . The inclusive reaction (two- and 
three-fragment) has been studied recently by Efros et al. 
p6| with realistic two-body interactions and three-body 
forces as input in the energy range up to 100 MeV by 



employing the Lorentz integral transformation method. 
Other recent theoretical work was devoted to polarization 
observables for p-d capture [^7|-^l| using realistic inter- 
actions. A discussion of polarization observables will be 
published in a subsequent paper. 

In Rcfs. p2|-p^|, we have treated the 3 He and 3 H pho- 
todisintegration and the inverse processes, i.e., the radia- 
tive capture of protons or neutrons by deuterons, within 
the integral equation approach discussed below. These 
calculations were based on the Bonn A, Bonn B and 
Paris potentials in Ernst-Shakin-Thaler (EST) expan- 
sion: Bonn A (EST), Bonn B (EST), and Paris (EST) 
p5| , p6| . We have demonstrated, in particular, the role of 
E2 contributions, meson exchange currents, and higher 
partial waves. A noticeable potential dependence was 
found in the peak region, i.e., for E x < 20 MeV. But 
it was also shown that the different peak heights are 
strongly correlated with the different three-nucleon bind- 
ing energies obtained for the potentials employed. The 
possibility of using the magnitudes of the cross sections 
as an independent test of the quality of the potentials, 
thus, appears rather restricted. 

The aim of the present paper is to extend the investi- 
gations of Refs. to photon energies from 20 MeV 
to 100 MeV. Aside from certain energy points there are 
up to now no other theoretical calculations for the differ- 
ential and the total cross section available in this energy 
range using the the full final-state amplitudes, realistic 
interactions and taking into account El and E2 contri- 
butions of the electromagnetic interaction. 

This paper is organized as follows: In section || we 
present the theoretical framework of our calculations. 
The results are discussed is section [II. Our conclusions 



are summarized in section [V 



II. THEORY 

The amplitude for the two-fragment photodisintegra- 
tion of 3 H or 3 Hc into a deuteron ijj^ and a neutron or 
proton of relative momentum q is given by 
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M(q)= ^i(q;Vrf|^cm|*Bs)5 

- ( -i(*IJJ era l* BS ) s . (1) 

Here I^E'bs)^ represents the incoming three-nucleon 

bound state while ^ J (q; | denotes the hnal contin- 
uum (scattering) state with outgoing boundary condi- 
tion. H cm is the electromagnetic operator. The initial 
and final states are assumed to be properly antisym- 
metrized. The antisymmetrizcd final state can be rep- 
resented as a sum over the three possible two-fragment 
partitions (3 

C 1(q;^|= ( 1(*I = ^E (-) <*/*l- ( 2 ) 

The scattering state ^{^p\ can be obtained from the 
free channel state (3>p\ via the M0ller operators fil \ 
i.e., ^(^p\ = {®p\ fi^*. The channel state 

= ((p)qSM s ;TM T \ (((3)Wmf,mt\ (3) 

is a tensor product of the deuteron wave function and a 
plane wave state of the outgoing third particle. More- 
over, it is an eigenstate of the channel Hamiltonian 
Hp = Hq + Vp. Here we have used the complementary 
notation Vp = V ia for the two-body potentials, i.e., V ia 
denotes the interaction between the particle 7 and a, 
while Hq denotes the free three-body Hamiltonian. The 
states in Eq. (^) are labeled by the corresponding quan- 
tum numbers. The collective index rj — (sj; t) stands for 
the spin s, the total angular momentum j, and the isospin 
t of the two-body subsystem. The indices S and r denote 
the spin and isospin of the third particle, respectively. 

It can be shown ||47| that the adjoint M0ller operators 
satisfy the relation 



" = 5p a + Up Q {Ep + iO) G a {Ep + iO) , (4) 

where Up a are the usual AGS |H| operators, and G a is 
the resolvent of the channel Hamiltonian H a = Hq + V a . 



partial wave Bonn A (EST) Bonn B (EST) Paris (EST) 

5 5 5 

3 Si - 3 Di 6 6 6 

1 L> 2 4 4 5 

3 L> 2 4 4 5 

1 P 1 4 4 5 

3 Pi 4 4 5 

3 P 4 4 5 

3 P 2 - 3 F 2 5 5 7 



TABLE I. Ranks of the two-body partial waves of the 
Paris, Bonn A, and Bonn B potentials in EST representa- 
tion used for the bound-state and the scattering calcula- 
tions. 



Multiplying the AGS equations 

Up a = (1 - SpajGy 1 + 2J(1 - 5p a )T 7 G U ia (5) 
7 

from the right with G a and adding 5p a on both sides, we 
see that the left-hand side is already the expression (0), 
and inserting this relation on the right-hand side of (|5|) 
we end up, after some trivial manipulations, with the set 
of integral equations 

^=1 + ^(1-^7)^7^00^ (6) 
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for the adjoint M0ller operators. Applying the M0ller op- 
erator (||) onto the state -Heml^Bs), they go over into a 
set of effective two-body equations when representing the 
input two-body T operator in separable form. In order to 
accomplish this, we use the separable expansion method 
proposed by Ernst, Shakin, and Thaler [4§| ] for represent- 
ing a given NN interaction. In this scheme the original 
potential is expressed as sum over separable terms 

JV 

vfl,= \9 v J)^Ag2i'\, (7) 
1 

where TV is the rank of the separable representation, A M „ 
are the coupling strengths, and \g% I) are the form factors. 
Here I and I' are the orbital angular momenta. The total 
angular momentum j is obtained from the coupling se- 
quence {ls)j. Using this representation for the potential, 
the two-body T operator reads 

Tl{E + *0) = ]T \gl I) A^(E + »0) (g2 1'\ , (8) 

ft is 

with 

A"(£ + zO)= ((A 1 !)- 1 - g (E + iO))' 1 (9) 

and 

(Go(E + iO))^ =J2(9n l \Go(E + iO)\g„l). (10) 
l 

For more details of this construction we refer to Refs. 
P5| , p6[ . The ranks for each partial wave used in this 
paper for the bound-state and the scattering calculations 
are contained in Tab. [j]. The separable representation 
reproduces the correct negative-energy bound-state pole 
of the two-body T matrix (for the deuteron quantum 
numbers r]d), if the form factor satisfies the homogeneous 
equation 

*Ev3,'G {E d )\g? l') = \g?l). (11) 
v 

Then the deuteron wave function is given by 

\i; d )='£G (E d )\gri). (12) 
1 
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It should be pointed out, that we have renormalized the 
form factors and coupling strengths in the 3 Si- 3 -Di chan- 
nel from the original representation in Eq. (]?]) in order to 
give the correct normalized deuteron wave function ac- 
cording to Eq. (|l2|) . This procedure does not change the 
potential or the T matrix, since we shift only a normal- 
ization constant from the form factors into the coupling 
strengths. 

Equation (^) will be treated numerically in momentum 
space, employing a complete set of partial wave states 
\pqlbTM r ;IMi). The label b denotes the set (rjSKL) 
of quantum numbers, where K and L are the channel 
spin of the three nucleons [with the coupling sequence 
{jS)K\ and the relative angular momentum between the 
two-body subsystem and the third particle, respectively, 
r is the total angular momentum following from the cou- 
pling sequence (KL)T, the total isospin / follows from 
the coupling (tr)I. These states satisfy the complete- 
ness relation 



i = EEE d*?** 

TM T lb IMi q q 

x \pqlb TM r ; I Mi) (pqlb TM r ; I Mi \ 



(13) 



The required antisymmetry under permutation of two 
particles in the subsystem can be achieved by choosing 
only those states which satisfy the condition (— )'+ s +* = 
-1 . 

Using the above defined states, the partial-wave de- 
composition of the channel state defined in Eq. (||) reads 

(^l = EE E Y. Y ^ti){jrn ] SM s \KM K ) 

TM r b M K M L IM, 

x {KM K LM L \TM T ){tmtrm T \IMi) 

x ((p) gi qbTM r ;IMi\G (E d +lq 2 + iO), (14) 

with 

((fi) gi qbTMr;IMi\ 

OO 

= E / d PP 2 (((3)pqlbTM T ;IMi\ g$(p) , (15) 



where we have used Eq. (|12|) for the representation of 
the deuteron wave function. As the quantization axis we 
have chosen the direction of the incoming photon. The 
generalization of Eq. ( |l5| ) for arbitrary rank index fi and 
arbitrary channel quantum numbers is given by 

m g)1 qbTMr;IMi\ 

OO 

= J2 [d PP 2 mpqlbTMr;IMi\ $(p). (16) 



These states are needed to derive and solve an integral 
equation for the final-state amplitudes. From the solu- 
tion of this equation only those channels with a deuteron 



as subsystem, i.e., the channels from Eq. (|l5|), are 
needed for the calculation of observables for the two-body 
breakup or the capture process. 

Equipped with the above equations we can now derive 
an integral equation for the scattering amplitudes. Mul- 
tiplying Eq. @ from the left with Gq and the states of 
Eq. (|l6|), and from the right with .Heml^Bs)^ we obtain 



ri AMl{ q ,E d + \ q 2 ) 



^AB h Jq,E d + h 2 ) 



(17) 



EE AW 2 TI AV^M,E d +\q 2 ) 
b' vp i 



xA? p (E d + lq 2 -lq'Y I AM b ;(q',E d 
with 

ri AMl(q,E) 



(18) 



4= E^) ^ qbr- I\G Q (E + iO)n<j )f H eia \* BS ) s , 
V3 p 



and 
ri AK! b 



<ABl(q,E) (19) 

= ]T<(/3) 9^Q9bT; I\G (E + i 0)iJ cm |tf BS )s , 
V3 

where we have used the separable expansion of Eq. (|^) 
for the T operator. Here TI AB b ll represents the so-called 
plane- wave (Born) amplitude and rl AM b /1 denotes the full 

final-state amplitude. The effective potential TI AV hb v en- 
tering Eq. ( fi"7| ) is given by 

ri AV%(q,E)=J2(l-S ,) (20) 



x (09) g, qbT; I\G (E + i0)\ g v q' b' V; J'( 7 )) . 

Noting that the two non-zero contributions to the sum 
are identical, we can write the effective potential as 



l AV%{q,4,E) = 2^ / / dpp 2 dp'p' 2 gl(p) (21) 



"00 

x mpqlbT; I\G (E + iO)\p' q' l' b' T; I (a)) <#'>')■ 

The recoupling coefficients entering this equation can be 
found in Ref. fdll (or in a more compact form in p^j for 
another coupling sequence which can easily be changed to 
the present one). In Eqs. (17)-(|l]) we have used the fact 
that the Born term, the effective potential, and therefore 
the full amplitude are diagonal in the quantum numbers 
T, M r , I, and Mi. 

In order to be able to solve Eq. (0) numerically an off- 
shell extension is required. This can easily be achieved 
by replacing E d + \q 2 with the energy parameter E. The 
solution of Eq. (p~7]) is obtained on the real axis by ex- 
panding the solution in cubic B splines pa and solving 
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a system of linear equations for the unknown coefficients. 
The logarithmic singularities in the kernel of the integral 
equation have been treated by a standard subtraction 
technique. For more details on the numerical solution of 
an integral equation of similar type we refer to Jl9| and 
references therein. 

In the Born amplitude one finds that the terms in the 
summation are independent from the partition, i.e., the 
summation over the different clusters can be replaced by 
a factor of 3. Using the states of Eq. ( |l6| ) we obtain 



ri ABl(q,E) = V3Yl J d PP 2 



E -p 2 - |<7 2 + i0 

xipqlbT-JlH^B^s. (22) 

For energies E above the deuteron breakup we have to 
take care of the pole in the propagator, which is done by 
using a standard subtraction technique. 

The three- nucleon bound state |\&bs) with binding en- 
ergy Ebs which is contained in the expression for the 
Born amplitude is determined by the Schrodinger equa- 
tion 

{Ebs - H) |*bs> = 0, (23) 
where the total Hamiltonian H is given by H = Hq + V = 

3 

Hq + Vy. When introducing the channel resolvents 

7=1 

Gy(z) = (z — H Q — V^) -1 , Eq. ( |23| ) can be written in 
form of a homogeneous integral equation, 



|* BS ) =G 7 (£ BS )^ 7 |1'bs> 

= G 7 (£bs)V(1-M^I*bs>, 



(24) 



with V~ ( — V — Vy being the channel interaction between 
particle 7 and the (a(3) subsystem. If we now introduce 
the position \Fp) = (V — Vp) I^bs), and use the relation 
VyG~/ = TjGq, we obtain the equation 

\F P ) = 2(1 - ^ 7 )T t (^bs)Go(^bs) \F- ( ), (25) 



where Go again is the resolvent of the free Hamiltonian. 
In Eq. (p5j), the summation runs over all two-fragment 
partitions 7. The "form-factors" \Fp) are related to 
|*bs> by 

I* bs) = 2 G o(£bs) T 7 (£ B s) G (Ebs) \FJ 

l^7> . (26) 



7 

7 



Bonn A (EST) 


Bonn B (EST) 


Paris (EST) 


-8.284 


-8.088 


-7.3688 



TABLE II. Calculated binding energies for the three- 
nucleon bound state. The total angular momentum of the 
two-body subsystem was restricted to j < 2. 



where the |^ 7 ) are the standard Faddeev components. 
For a numerical treatment of Eq. ( |25| ) we multiply this 
equation with Go and the partial- wave states (p q I b T; I\ . 
After inserting the separable T matrix from Eq. (ph and 
defining 



(?) = E/ dpp 2 9l(j>)(pqlbT-I\G o \F ), (27) 



Eq. 



goes over into 



^W-EE jdq'q' 2 AV%(q,q\E BS ) (28) 
b' vp { 

xA^EBS-^F^iq'). 



After discretization Eq. (|28| ) can be treated as a lin- 
ear eigenvalue problem, where the energy is considered 
as a parameter which is varied until the corresponding 
eigenvalue equals unity. The eigenvalues can be found 
by using standard numerical algorithms. The obtained 
binding energies for the three potentials used in this pa- 
per can be found in Tab. 0. As shown in 52 these 
values are practically the same as those for the original 
potentials. 

The whole wave function can now be calculated by 
either using Eq. (|2^), or by applying the permutation 
operator P on one Faddeev component |5Q] 



|*BS) = (1 + P)|^l), 



(29) 



where P represents the sum of all cyclical and anticyclical 
permutations of the nucleons. The required antisymme- 
try of the wave function is achieved by projecting only 
on those channels with (— )'+ s +* = — 1. 

In our calculation of the Faddeev components the to- 
tal angular momentum j of the two-body potential was 
restricted to j < 2 (18 channels), while in the full state 
all partial waves with j < 4 (34 channels) have been 
taken into account. With this number of channels con- 
verged calculations of the observables for the photopro- 
cesses in this paper were achieved, incorporating 99.8% 
of the wave functions j|2],fl3| . For more details concern- 
ing the properties of the wave functions and their high 
quality we refer to |5g] . 

The relevant electromagnetic operator in the total 
cross section at low energies is a dipole operator. In the 
differential cross section and at higher energies also the 
quadrupole operator is relevant. In order to take into ac- 
count meson exchange currents we use Siegert's theorem 
p3[ , then these operators are given by fl54fl 



H 0) 




(30) 
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FIG. 1. Total cross section for the photodisintegration of 
3 He. The data are from [2-4,8-10,12-14,20]. 
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FIG. 3. Differential cross section at 90° for the photodisin- 
tegration of He. The data are from [2,5,8-13]. 



and 




eir?y 2 A(^, <Pi), 



(31) 



where E 7 denotes the photon energy. Here, rj are the 
nucleon coordinates, the electric charges, and A = ±1 
is the polarization of the photon. The normalization fac- 
tor M contains the quantization volume and is canceled 
out in the calculation of the cross sections. For the cal- 
culation of matrix elements of these operators with ini- 
tial and final state they need to be transformed into the 
three-body center-of-mass system. Expressions for the 
matrix elements can be found in ]55| ]. 

The on-shell amplitudes can be obtained from 



(-) 



"j < q S M s ; ipdjmj | H cm | * BS r'Mrv ) s (32) 

2 

EEE E Y LML {q)UmjSM s \KM K ) 

r =i M T b M K M L 



x ( KM K LM L | TM r ) TI AM\ (q, E, 



h 2 ) 



In the summation of course only those channels con- 
tribute that have a deuteron as their subsystem. The 
total angular momentum of the bound-state wave func- 
tions is r = i. This means when taking into account 
E2 contributions, the maximum total angular momen- 
tum of the outgoing state can be T = §. With these 
amplitudes the unpolarized differential cross section for 
the photodisintegration process is given by 



x E E | ( ~i ( <1 S Ms > I H ^ I * BS r ' Mr ' 

M s M.j AM' 



(33) 

2 



The differential cross section is usually expanded in terms 
of Legendre polynomials 

<r(q, 0) - ^( ? , 6) = A (l + £j a k P,(cosr5)^ . (34) 

The coefficients A and can be calculated analytically 
from Eqs. ([32]) and (|33|). The total cross section is ob- 
tained by integrating Eq. (B4) over the angle 9 between 
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FIG. 2. Total cross section for the photodisintegration of 
3 He. The data are from [4,9,10,16]. 
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H. The data are from [22-28]. 
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E x [MeV] An [fib] ai a :i a4 



91 47 


U.O i o 


n 49 


-u.oo 


-U.4J. 


-u.uo 












fl D74-D 1 


OA 1 A 


U.ol ( 


C\ A A 
U.44 


-U. (\) 


n a q 


-U.Uo 




0.42± 0.04 


0.34±0.1 


-0.86±0.12 


-0.35±0.15 


-0.12±0.1 


26.81 


0.467 


0.46 


-0.75 


-0.44 


-0.06 




0.38± 0.03 


0.41±0.1 


-0.86±0.13 


-0.39±0.15 


-0.11±0.1 


29.47 


0.425 


0.47 


-0.72 


-0.45 


-0.06 




0.33± 0.03 


0.39±0.1 


-0.77±0.10 


-0.42±0.15 


-0.12±0.1 


32.14 


0.389 


0.48 


-0.68 


-0.45 


-0.07 




0.33± 0.03 


0.33±0.1 


-0.79±0.12 


-0.36±0.15 


-0.08±0.1 



TABLE III. Coefficients for the expansion of the differenti 
corresponds to the theoretical results obtained with the Paris 
Ref. [19]. 

the incoming photon and the outgoing proton or neutron 



= 47T An 



(35) 



The cross section for the p-d or n-d capture process 
is obtained from the corresponding photodisintegration 
expression by using the principle of detailed balance |35|] 



do 



(.lis 



dn 



3 k 2 der ca P 
2 Q2 



dn 



(36) 



Here, k and Q are the momenta of the proton and 
the photon, respectively. In the present treatment no 
Coulomb forces have been taken into account, in other 
words the matrix elements of Eq. ([!]) for p-d capture dif- 
fer from the corresponding n-d expression only in their 
isospin content. 



III. RESULTS 

It should be pointed out that we have shifted all the- 
oretical cross sections to the experimental threshold for 
a meaningful comparison. All calculations are done with 
the theoretical binding energies. 
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Total cross section for the capture of protons by 
The data are from [1,6-8,12,13,18,19,21]. 



al cross section for p-d capture. For each energy the first row 
(EST) potential, and the second row contains the data from 

In Fig. [l] we show our theoretical results for the to- 
tal cross section of 3 He photodisintegration compared to 
most of the available experimental data P-|^J^-[l3]] up to 
E 1 = 40 MeV. Not shown are the data by van der Woude 
et al. |l5| and Chang et al. JlTj who found evidence for an 
excited state in their measurements. This resonance be- 
havior has never been confirmed by any other group. It 
can be seen from Fig. [ljthat there are large discrepancies 
between the different data sets around E 1 — 1 1 MeV. The 
theoretical curves lie in between the data sets. It should 
be emphasized that for the calculated curves there is a 
correlation between the three-body binding energy and 
the peak height of the cross section for the photodisin- 
tegration p4p3[ , i.e., the higher the binding energy the 
lower the cross section at the peak. Above 12 MeV the 
mentioned discrepancy of the experimental data declines. 
Due to the large error bars it is not possible to draw fur- 
ther conclusions. 

In Fig. H we present total cross section calculations of 
photodisintegration of 3 Hc between E~ f = 40 MeV and 
E 1 — 100 MeV compared to the measurements of Feti- 
sov et al. [Q, Kundu et al. Q, Ticcioni et al. ]Io[ ), and 
O 'Fallon et al. @. For energies above E^ = 60 MeV 
the measured points lie slightly above our curves, com- 
puted by employing Bonn A, Bonn B, and Paris poten- 
tials. However, our curves agree with the tendency of the 
data. 

We would like to point out that especially in this en- 
ergy range a high rank representation of the NN poten- 



E x [MeV] 


A [fib] 


ai 


ai 


az 


a 4 


21.47 


0.579 


-0.071 


-0.88 


0.069 


-0.0021 


24.14 


0.517 


-0.065 


-0.84 


0.064 


-0.0023 


26.81 


0.464 


-0.057 


-0.81 


0.056 


-0.0026 


29.47 


0.421 


-0.047 
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-0.035 


-0.75 


0.037 


-0.0030 



TABLE IV. Coefficients for the expansion of the dif- 
ferential cross section for n-d capture for the Paris (EST) 
potential. 
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FIG. 6. Angular distribution and differential cross section for the capture of protons by deuterons. The data are from 



[19,21]. 

tials is required in order to get converged results. Below 
E 1 = 40 MeV the improvements with respect to a low 
rank calculation are of the order of 1-5%. In view of 
the experimental errorbars this change is of course not 
relevant. Above = 40 MeV the low rank calculations 
yield a cross section which is 5-15% lower than the high 
rank calculations presented in this paper. 

The differential cross section calculations at 90° for the 
photodisintegration of 3 He up to an energy of E 7 = 40 
MeV are illustrated in Fig. || in comparison to the corre- 
sponding experimental data ||J|J|-(l|] . The EST repre- 
sentations of the potentials Bonn A, Bonn B, and Paris 



are employed. As in the previous figures, we can not ob- 
serve any significant difference for energies above E 1 = 20 
MeV, whereas the peak region shows a considerable po- 
tential dependence, as discussed for Fig. [l]. The data by 
Berman et al. j^] are below the calculated curves, how- 
ever, they do agree with the tendency in the peak region. 
There is a remarkable discrepancy between these data 
sets and the data points measured by Kundu et al. [||, 
but they coincide with the theoretical curves for energies 
above 25 MeV. We find again that the theoretical results 
lie in between the data, though there are discrepancies 
between the data sets and, moreover, the errorbars are 
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FIG. 7. Angular distribution and differential cross section for the photodisintegration 3 He. The data are from [16] 



quite large. 

The two-body photodisintegration of 3 H has been mea- 
sured by Bosch et al. ]22]], Kosiek et al. |23|p4|], Faul et 
al. and Skopik et al. j2(|. In Fig. |we display 

these data compared to our theoretical calculations. Also 
shown in this figure are the transformed results by Mitev 
et al. (27]] and Mosner et al. |||| obtained from radia- 
tive n-d capture measurements. It should be pointed out 
that the most recent measurement by Mosner et al. is 
in excellent agreement with the theory. We notice that 
for low energies, e.g., up to E 1 — 15 MeV, the calcu- 
lated curves for the different potentials show a different 
behavior, whereas for higher energies all three calcula- 
tions do not yield any significant difference. In the peak 
region, the curves for the Bonn A and the Bonn B po- 
tentials cover the experimental data in between the error 
bars better. For energies above 20 MeV there is a large 
discrepancy between the data sets by Kosiek et al. and 
Skopik et al., although the tendency of the data is similar. 
This indicates a normalization problem. 

In Ref. jlH we discussed the available data for p-d cap- 
ture below E x = 20 MeV. It was shown, that only the 
coefficient Aq of the expansion in Eq. (|34|) has some po- 
tential dependence, whereas the coefficients are almost 
independent from the interaction. Also in this case there 
is a correlation between the peak height and the binding 



energy, i.e., the lower the binding energy the lower the 
peak height. It should be pointed out, that this is the 
inverse of the relation found in case of the photodisin- 
tegration. We also have demonstrated that there seems 
to be a normalization problem in the experimental data. 
It can be seen in Fig. @ that the data by Wolfli et al. 
IIkJ, Matthews et al. [|l8| and Anghinolfi et al. |lj| are 
too low compared to those by King et al. M], Pitts et 
al. [f21[ , and Belt et al. || which agree with our theo- 
retical curves. This indicates a calibration problem of 
the measurements. It was also shown in [Q that after 
renormalization the data sets by Matthews et al. are in 
agreement with those by King et al. and the theoretical 
curves. At energies above E 1 — 20 MeV we encounter 
a similar problem and compare in Fig. ^| the differen- 
tial cross section divided by A$. It can be seen that the 
agreement between theory and the experimental data by 
Anghinolfi et al. is very good. A comparison of the ex- 
pansion coefficients obtained by Anghinolfi et al. and our 
theoretical values for the Paris (EST) potential is given 
in Tab. 



[II 



There are discrepancies for Aq which are 
connected to the normalization problem mentioned ear- 
lier. Despite the relatively big experimental error bars 
for the expansion coefficients there are considerable 
good agreements. 

Also shown in Fig. || are the data by Pitts et al. |plj . 
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FIG. 

by deuterons at E^ h 



In this case there is also excellent agreement for the ab- 
solute cross section, particularly by employing the Bonn 
A potential. There are two additional data sets by van 
der Woude et al. at E x = 19.2 MeV and E x = 20.6 
MeV which are not shown here because of the measured 
resonance behavior, as mentioned above. 

Besides total cross section data, O 'Fallon et al. 
have also measured data sets of the differential cross sec- 
tion for photodisintegration of 3 He up to an energy of E 1 
— 140 MeV. As can be seen in Fig. || their total cross sec- 
tions are slightly higher than the theoretical predictions. 
For a meaningful comparison with our calculations we il- 
lustrate in Fig. [?] four of their data sets for the differential 
cross section normalized with Aq. Within their error bars 
they agree quite well with the theoretical calculations. 

The different contributions for the El and E2 tran- 
sitions in case of p-d capture at E^ ab = 10.93 MeV are 
shown in Fig. ||. The pure E2 contributions are very 
small and enter the differential cross section essentially 
through the E1-E2 interference term, which leads to the 
asymmetry, i.e., the curve is shifted to smaller angles. 
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FIG. 9. Differential cross section for the capture of neu- 
trons by deuterons at E^ h = 10.8 MeV for the Bonn B po- 
tential. The data are from [27]. 
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FIG. 10. Differential cross section for the capture of neu- 
trons by deuterons. The data are from [27]. 

With the inclusion of E2 transitions there is an excel- 
lent agreement with the data by King et al. fll3fl . The 
only measurement of the angular distribution of the dif- 
ferential cross section for n-d capture has been done by 
Mitev et al. (27]]. The different contributions for the El 
and E2 terms is shown in Fig. || for E l * b = 10.8 MeV. 
Due to isospin selection rules the E2 contribution, and 
hence the interference between the El and the E2 term, 
is much smaller compared to 3 He. It should be noted that 
in this case the maximum of the differential cross sec- 
tion is shifted to larger angles. This observation was also 
made in Ref. |29| for the Born approximation. A compar- 
ison to the theoretical calculations at E\f° — 9 MeV and 
E\f h = 14 MeV is shown in Fig. [l0| It can be seen that 
the peaks of the experimental data tend to have a bigger 
asymmetry than the theoretical curves. It is remarkable 
that this circumstance is independent of the potential 
choice. This indicates either a stronger contribution of 
a higher multipole, not present in our theoretical cal- 
culations, or an error in the data. To the best of our 
knowledge there are no other differential cross section 
data for 3 H photodisintegration or the inverse reaction 
available. There are also no experimental data available 
for the Legendre coefficients dfe. Nevertheless, we show 
for comparison in Tab. IV corresponding calculated val- 



ues for n-d capture at the same energies as the available 
data for p-d capture from Tab. [II. It can be seen that 



for n-d capture the angular distribution is dominated by 
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text. The data are from [6-9,12,13,18-20]. 

the a-i coefficient, i.e., the El transitions, whereas in case 
of p-d capture there are bigger admixtures of E2 contri- 
butions. 

One observable of particular interest in the context of 
angular distributions is the so-called fore-aft asymmetry. 
This quantity is defined by 



cr(54.7°) - cr(125.3°) 



(37) 



cr(54.7°) +cr(125.3°) ' 

In terms of Legendre polynomial expansion coefficients 
of Eq. (34), this can be written as 



a i 



..... (38) 
V3(l - yga 4 ) 

In Figs, [ll] and [l2] we compare our results with the avail- 
able data sets. The theoretical curves for 3 He agree quite 
well with the data | |,^|jl|,|l|-|§ , whereas the calcu- 
lated asymmetry for H is smaller by a factor 5 than 
the experimental data |p2| , p6 27 j. A similar observation 
was made by Skopik et al. ]20|] using their effective cap- 
ture calculations, where no FSI effects were taken into 
account. Since all experimental data show a consistently 
higher fore-aft asymmetry, there seems to be something 
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missing in the theoretical description of this reaction. 
One possible explanation for the discrepancy between 
theory and experiment was that the FSI has not been 
taken into account properly in previous calculations. In 
the present calculations we have shown that this discrep- 
ancy still remains when taking FSI effects into account. 
Therefore, it is still unclear where the differences stem 
from. The Ml term is not likely to solve the problem, 
since it is only expected to have an effect at extreme an- 
gles or at very low energies. Also three-nucleon forces 
are not expected to solve the problem since the angular 
distribution shows no potential dependence. A possible 
solution could be the inclusion of explicit meson exchange 
currents which allow a stronger coupling of higher mul- 
tipolcs to 3 H. 

IV. CONCLUSIONS 

In this paper we have analyzed all available experimen- 
tal data for the photodisintegration of 3 He and 3 H and 
the corresponding inverse reactions below E 1 = 100 MeV 
by comparing with our calculations using realistic NN 
interactions. We have shown that the theoretical curves 
agree with the experimental data for the total cross sec- 
tion within the error bars. Moreover, in many cases the 
measured differential cross sections for p-d capture (aside 
from a normalization factor) can be explained theoreti- 
cally over a large energy range. In Q] it was already 
shown that a similar normalization problem exists for 
the data below E x = 20 MeV. There, it was also shown 
that the angular distribution is insensitive to the under- 
lying two-body interaction, whereas there is a strong cor- 
relation between the three-body binding energy and the 
normalization constant Aq [ fl2|]43]] . Since the angular dis- 
tribution is insensitive to the employed interaction, we 
do not expect large effects of three-nucleon forces. On 
the other hand taking account of them will change the 
three-body binding energy and hence the normalization 
constant Aq. 

For n-d capture the description of the angular distri- 
bution is less good. For energies above 10 MeV the the- 
oretical results give a much smaller asymmetry than the 
experimental data. Hence, the theoretical fore-aft asym- 
metry shows a large discrepancy from the experimental 
data, whereas for p-d we achieve a very good agreement. 
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